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INTRODUCTION 


Thermomechanical service conditions for high-temperature structural components, 
e.g., hot-section aircraft engine and nuclear reactor components, involve tempera- 
ture levels, thermal transients and mechanical loads severe enough to cause measur- 
able inelastic deformation. Although inelastic strain per cycle is small in many 
applications (in the order of elastic strain) as cycling occurs in the absence 
of shakedown inelastic strain can accumulate monotoni cal ly (ratchetting) or alter- 
nate repeatedly (fatigue). In any case, the seat of failure clearly resides in 
the occurence of continued inelastic deformation. Structural analysis in support 
of the design of high-temperature components - leading to the stress, strain and 
temperature fields upon which life predictions are ultimately based - therefore 
depends strongly on accurate mathematical representations (constitutive equations) 
of the nonlinear, hereditary, inelastic behavior of structural alloys at high 
temperature, particularly in the relatively small strain range. To be generally 
applicable, constitutive equations must be expressed in multi axial form and be 
appropriate for all modes of mechanical and thermal loading to be experienced by the 
elevated temperature components (e.g., cyclic, non-isothermal , non-radial, etc.). 

In contribution to the overall program in constitutive equation development 
sponsored by the HOST project, the in-house NASA/LeRC effort is concentrating mainly 
on three fundamental areas of research: 

1) Multiaxial experimentation to provide a much needed rational basis for high- 
temperature multiaxial constitutive relationships. 

2) Noni sothermal testing and theoretical development toward a complete thermo- 
mechanically path dependent formulation of viscoplasticity - including the 
influence of thermally induced metallurgical changes. 

3) Development of viscoplastic constitutive models accounting for strong ini- 
tial anisotropy. 

Current activities, progress and future directions in these areas of research 
are indicated in the following sections. 


* Work performed under NASA Grant NAG-3-379. 
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MULT I AX I AL EXPERIMENTATION IN SUPPORT OF CONSTITUTIVE EQUATION DEVELOPMENT 

Thoughtful exploratory experimentation must go hand-in-hand with the formula- 
tion of constitutive theories and furnish guidance for their development. Without 
close interaction between experimentalist and theoretician physically unrealistic 
and ad hoc constitutive models often result that may not be used with confidence 
even slightly outside the specific conditions addressed by their data base. 

A particularly glaring deficiency in terms of elevated temperature testing is 
that regarding multiaxial behavior. A major reason for this has been the lack of 
extensometry for accurately measuring multiaxial strain at high-temperature. Only 
through recent developments in extensometry (e.g., ref.l) has it become possible to 
do meaningful multiaxial testing at relevant temperatures on relevant materials. 

The necessity of conducting fundamental multiaxial tests is best recognized in 
the context of classical plasticity theory. In that constitutive theory the con- 
cept of a yield surface plays a central role. A description of the yield surface, 
at a fixed inelastic state, together with the concept of normality is precisely 
that which allows a consistent multiaxial plastic flow law to be expressed. An 
appropriate description of the yield surface and a demonstration of the validity of 
the concept of normality can only be realized through multiaxial testing. 

At high temperature, alloys of interest are strongly time-dependent (visco- 
plastic) and the concept of a yield surface, in the classical sense, breaks down. 
However, concepts have been postulated for time-dependent inelastic behavior that 
have an analogous geometrical interpretation to that of yield surfaces (e.g., 
Odqvist (ref. 2) and Drucker (ref. 3)). Such a concept was introduced by Rice (ref. 
4), Ponter and Leckie (ref. 5) and Ponter (ref. 6) in the form of a flow potential 
function of the applied stress o-jj and the internal (back) 

ii (°i j > a i j ) 

from which the viscoplastic flow and evolutionary laws are 

• _ d\l 

£ ij ■ ao-jj 

«i j = -h(«kl ) 3ai j 

Here, e-jj is the inelastic strain-rate and h is a scalar function of the internal 
stress. Experimental evidence for the existence of a potential function was found 
in a preliminary way by Brown (ref. 7) and Robinson (ref. 8) and recently and more 
comprehensively by Oytana, Delobelle and Mermet (ref. 9). 

Robinson (refs. 10,11 ) has adopted a potential function of the general form:** 

u = 2 Ijf(F)dF + jj-/g(G)dG (4) 

* As the tests described here are isothermal, the development in this section is 
similarly limited to isothermal conditions. Also, isotropic hardening effects are 
taken to be saturated out. 

** The complete model utilizes a discontinuous pontential function that allows 
analytically different governing equations in various "regions of the state space. 
These aspects of the model will not be discussed here. 


stress aij 


derivable, i.e. 


( 1 ) 


( 2 ) 


(3) 
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in which o, R and H are material constants (or, under non-i sothermal conditions, 
functions of temperature) and the stress dependence enters through the functions 

F (°i j » u i j ) and G(ai j ) (5) 

The flow equation (2) then becomes 

2 “ £ u = f < F >4" < 6 > 

which can be interpreted geometrically in a superimposed stress and inelastic 
strain-rate space (fig. 1) as indicating normality of the inelastic strain-rate 
(vector) to the hypersurfaces 

F (°i j , oti j ) = const. (7) 

The internal stress a-jj, denoting the inelastic state, is taken fixed. 

The surface F = 0 is considered by Robinson to correspond to a threshold 
(Bi ngham-Prager) stress below which the inelastic strain-rate vanishes. 

Experimental determination of the function F is analogous to the determination 
of the yield function (or a yield surface) in classical plasticity. Specification 
of F at constant inelastic states permits the definition of a consistent multi- 
axial form of the viscoplastic flow law, and provides information concerning an 
evolutionary law for the inelastic state variable a-jj. 

For isotropic alloys it is reasonable to further specialize F and G to be 
dependent on the principal invariants of the applied and internal stress. In the 
spirit of v.Mises, Robinson has taken 

J 2 - ^2 

F(J 2 ) = — 1 and G(J 2 ) = — — (8) 

K2 k2 

as depending on only the second principal invariants 

J 2 = 7 £ ij £ ij (9) 

J 2 = \ a ij a ij U°) 

The first of equations (8) plays the role of a (Bingham) yield condition with K 
denoting the magnitude of the threshold Bingham shear stress. 

In equation (9) the effective stress 

£*ij = s ij" a ij 

is taken as the difference of the applied deviatoric stress S-jj and the internal 
stress deviator a-jj. The choice of stress dependence in equation (8) is consistent 
with the experimental observations reported in reference 9. 
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With these choices. 

equations (2) and (3) give 
2u£ij = f(F)E-jj 

(12) 

and 

(13) 

• • 

a i j = h(G)e-j j-r(G)a*j j 

with 



F- ^ (G) 

(14) 

Evidently, the evolutati onary equation (13) is of the physically accepted 


Bai 1 ey-Orowan form. 

of the equation (12) we get 


Squaring both sides 



4vi2j = f2j 2 = p(J 2 ) 

(15) 

where 

T 1 * 

1 ' 2 e ij e ij 

(16) 


is a measure of the magnitude of the inelastic strain-rate. Equation (15) states 
(see Fig. 2) that surfaces of constant J 2 in stress space (at a fixed inelastic 
state) are surfaces of constant inelastic strain-rate I. Further, the strain-rate 
vectors are normal to these surfaces. Constant inelastic strain-rate surfaces are 
the counterparts of yield surfaces (surfaces of constant inelastic strain) in time- 
dependent plasticity. 

In general, these surfaces change position and size through changes in the 
internal stress, threshold stress and time as inelastic deformation and recovery 
occurs. A complete understanding of the nature and behavior of these surfaces 
under virgin conditions and conditions subsequent to inelastic deformation (creep, 
cyclic plasticity, etc.) is prerequisite to formulating a consistent description of 
multiaxial viscoplastic behavior. 

Computer controlled experiments for directly determining the loci of constant 
inelastic strain rate at fixed inelastic state are included in the test plan for 
the currently expanding Structures Division laboratory. In the interim period, 
before that laboratory if fully operational, preliminary tests of this type are 
being conducted under subcontract to Oak Ridge National Laboratory. 


NONISOTHERMAL TESTING AND THERMOMECHANICAL MODELING 

Essentially all of the important structural problems related to the design of 
aircraft engine hot-section components are noni sothermal . Inelasticity does not 
occur in these components as a result of mechanical loading alone but is generally 
induced through thermal transient cycles and thermal gradients. Nevertheless, the 
constitutive equations and life prediction models used in structural analysis and 
design are almost always based completely on experimental data collected under 
isothermal conditions. Isothermal tests are commonly conducted over the tempera- 
ture„and then "fit" as functions of temperature across the temperature range to 
furnish a "nonisothermal " representation. This traditional approach leads to "non- 
i sothermal" models that do not reflect the strong thermomechanical path dependence 
observed in, for example, the cyclic hardening behavior of some alloys of inter- 
est - particularly in the presence of metallurgical changes. 

Cyclic hardening of some common structural alloys within their temperature 
range of interest is believed to be influenced by the phenomenon of dynamic strain 
aging. Strain aging occurs in solid solutions where solute atoms (e.g., carbon, 
nitrogen, etc.) are particularly free to diffuse through the parent lattice. It is 
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energetically preferable for these solute atoms to occupy sites in the neighborhood 
of mobile dislocations where their presence immobilizes the dislocations or at 
least makes their movement difficult, thus causing strengthening. 

Isothermal cycling at temperatures where such metallurgical changes occur might 
therefore be expected to show abnormal hardening, i.e., higher hardening rates and 
greater saturation strengths than at temperatures both lower and higher. Macrosco- 
pic evidence of strain aging in three common alloys (i.e., Hastelloy X and types 
304 and 316 stainless steel) is shown in figures 3 through 5. In each case the 
hardening rate and the stress range at "saturati on", peak at an intermediate tem- 
perature in the range. This hardening peak is interpreted as a manifestation of 
dynamic strain aging. At lower temperatures the mobility of solute atoms is far 
less and strain aging cannot occur; at higher temperatures normal recovery 
precesses, e.g., climb of edge dislocations, take over. 

In the aging process described dislocations can, under some circumstances, 
break away from their solute atmospheres becoming mobile again. Although tempor- 
arily freed, dislocations can again be immobilized as solute atoms gradually 
diffuse back to them. As the thermally activated process of diffusion is involved 
and solute atoms are migrating to dislocations which themselves are moving under 
the applied stress, it is expected that the ensuing inelastic deformation (cyclic 
hardening in particular) has a complex dependence on thermomechani ca 1 history. 

Phenomenological evidence of thermomechanical path dependence under cyclic con- 
ditions is seen in the results of the simple nonisothermal tests reported in 
figures 3 and 4 (dotted curves). In these tests cycling is initiated at one tem- 
perature and after some cycling the temperature is changed and cycling resumed. 

Figure 3 shows the results of two nonisothermal tests on Hastelloy X cycled 
over a strain range at constant strain rate. In one, the specimen is cycled at 
800F for five cycles; the temperature is then changed to 1000F and cycling is 
continued to apparent saturation, at about one hundred cycles. In the second, this 
history is repeated up to thirty cycles where the specimen is then brought back to 
800F and cycling continued. Results of a similar test on type 304 stainless steel 
are shown in figure 4 (dotted line). 

The features of these test results that reflect thermomechanical history depen- 
dence are : 1) The change in strength (stress range) with temperature at a fixed 
number of cycles is always negative, i.e., an increase in temperature always pro- 
duces a decrease in strength and vice versa, contrary to the implication of the 
isothermal data; 2) The current strength, in particular the "saturation" strength, 
depends on the temperature-strain history. Evidently, the information contained in 
the isothermal data is not sufficient for a complete nonisothermal description of 
the cyclic deformation in the temperature range of interest. In fact, the data 
suggests that, with accompanying metallurgical changes, the materials retain a full 
memeory of their thermomechanical history to cyclic saturation. Similar observa- 
tions have been made on two-phase alloys (ref. 12) presumably arising from micro- 
structural changes associated with precipitation of the y' phase. 

Characterization of isotropic cyclic hardening effects requires that a consti- 
tutive model contain a scalar state variable (e.g., K) and a corresponding evolu- 
tionary law describing its rate of change with thermomechani cal history. A 
proposed form of evolutionary law appropriate for cyclic stressing (in the absence 
of thermal recovery) is 

K = F(P,T)P + G(P,T)T (17) 

i n which . 

p = J (e P .e P .) 1/2 (18) 

ij 
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is a measure of the accumulated plastic strain and T is the temperature. The 
current value of K (or equivalently the stress range) is determined from equa- 
tion (17) only if the thermomechanical path is known. 

Information for characterizing the function F can be obtained from ordinary 
isothermal tests, however, noni sothermal tests trust be conducted to supply informa- 
tion about G. Candidate tests for this purpose have been identified and will be 
conducted as the noni sothermal testing facilities become available in the 
Structures Division laboratory. 

Using preliminary information about the function G an appropriate evolutionary 
law having the form of equation (17) has been adopted for use with the Robinson 
viscoplastic model described in references 10 and 11. Predictions made by this 
noni sothermal model give excellent qualitative agreement with the experimental 
results of figures 3 through 5. 


A VISCOPLASTIC MODEL FOR TRANSVERSE ISOTROPY 

The need for increased efficiency in energy systems places greater demands on 
the high-temperature structural alloys used for system components. As higher oper- 
ating temperatures are sought, advanced materials are being developed to meet these 
demands. Good examples of such materials are the directionally solidified poly- 
crystalline alloys finding application in turbine blades of aircraft engines. The 
directional properties of these metals render them highly anisotropic relative to 
conventional alloys. This introduces additional complexity in understanding and 
mathematically representing their mechanical behavior over and above the already 
enormous complexities associated with elevated temperature. 

Here, we discuss a means of extending the isotropic viscoplastic model dis- 
cussed in the first section and fully described in references 10 and 11 to be 
applicable to materials with initial anisotropy (ref. 13) arising from directionally 
solidified grain growth. 

The direction of grain solidification at each point can be characteri zed by a 
field of unit vectors d-j(x|<). The mechanical behavior must then depend not only 
on the stress and deformation history at the point but also on the local 
preferential direction. It follows that equation 5 must be replaced by 

F ( Z i j , di d j ) and G(a-jj,d-jdj) (19) 

In the case of full isotropy F and G were taken as depending on the (second) 
principal invariants of E-jj and a^j (eg. ( 8 )). Here instead, we take F and G to 
depend on invariants that reflect the appropriate anisotropy. The theory of 
tensorial invariants (ref. 14) requires that for form-invariance under arbitrary 
rigid-body rotations F and G must be expressible in terms of the principal invar- 
iants of their respective tensorial arguments and invariants involving various 
products of these tensors. Here, we take the functions F and G as depending on 
subsets of these invariants^ We denote I i , I 2 and I 3 as the invariant arguments of 
the function F and I^, I 2 , T3 as those for G. 

That is, we take 

Fdi.i2i3) - 4~ )(, 2 - - 1 < 20 > 

K 2 |C K 2 J 
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As before, the flow and evolutionary equations are obtained by taking the 
appropriate derivatives as indicated in equations (2) and (3). This results in 

2 “' s ij - flFHqj * <-£ ~ * ViV 7 VkV ( “ij + d i d J ))) <23) 


for the flow law, and 

a. . = h (G ) e . . - r(G ) [a. . + ( 


l) ( dj d k a ki + d kdi a jk 

" 7 d p d k a kp (6 ij + d i d j ) ^ 


for the evolutionary law. 

As earlier, K denotes the threshold (Bingham) shear stress transverse to the 
preferential material direction. denotes the same for shear along the material 
direction. For K = K^, indicating no difference in shear strength across and along 
the direction d-j , equations (23) and (24) reduce to their respective isotropic 
counterparts eqs. (12) and (13). 

Equations 20 through 24 have been applied in several structural problems in- 
volving creep and cyclic plasticity. The predictions give good qualitative agree- 
ment with existing experimental results on the directionally solidified (DS) alloy 
MAR-M-247 as reported in reference 14. Quantitative comparisons are delayed until 
characterization tests on a DS alloy furnish specific values of the material para- 
meters in the viscoplastic model. 

Additional studies are continuing concerning the adoption of other sets of in- 
variants (integrity bases) reflecting higher degrees of anisotropy, e.g., ortho- 
t ropy and cubic symmetry. 
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Figure 1 Flow potential surfaces in stress space 

showing normality of inelastic strain rate vectors. 



Figure 2 For an isotropic material surfaces of constant 
inelastic strain rate are surfaces of J 2 =• constant. 
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Figure 3 lsothermal(solid)and nonisothermal(dotted) 
cyclic hardening curves for Hastelloy X. 
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Figure 4 lsothermal(solid)and nonisothermaKdotted) 

cyclic hardening curves for type 304 Stainless Steel. 
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Figure 5 Isothermal cyclic hardening curves for type 
316 Stainless Steel. 
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